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Localised states are found in many pattern forming systems. The aim of this paper is to investigate 
the occurrence of oscillatory localised states in two-dimensional Boussinesq magnetoconvection. 
Initially considering an idealised model, in which the vertical structure of the system has been 
simplified by a projection onto a small number of Fourier modes, we find that these states are 
restricted to the low ( regime (where C represents the ratio of the magnetic to thermal diffusivities) . 
These states always exhibit bist ability with another non-trivial solution branch, in other words they 
show no evidence of subcritical behaviour. This is due to the weak flux expulsion that is exhibited by 
these time-dependent solutions. Using the results of this parameter survey, we locate corresponding 
states in a fully-resolved two-dimensional system, although the mode of oscillation is more complex 
in this case. This is the first time that a localised oscillatory state, of this kind, has been found in 
a fully- resolved magnetoconvection simulation. 



I. INTRODUCTION 

Localised states are found in a wide range of phys- 
ical systems. Examples include the localised buckling 
patterns that are observed in the bending of beams [1 , 
cavity solitons in a semiconductor-based optical amplifier 
[2], isolated current filaments in a gas discharge system 
[3] and soliton-like structures on the surface of a ferrofluid 
[4]. Oscillatory localised states ('oscillons') have been ob- 
served experimentally in vertically- vibrated media |5l |6] . 
With appropriate choices for the nonlinear terms, the 
Swift-Hohenberg equation exhibits a range of localised 
stable solutions in both one [7-10 and two spatial di- 
mensions [12] , whilst equations of Ginzburg-Landau 
type can also produce localised states [13l |14] . Spatially- 
localised behaviour has also been found in simulations 
of two-dimensional thermosolutal convection p!5HT8] as 
well as in two-dimensional simulations of convection in a 
binary fluid [ME]. 

The presence of a strong vertical magnetic field 
tends to inhibit convective motions in an electrically- 
conducting fluid [22^ . This phenomenon can be observed 
on the surface of the Sun, where the dark central part of 
a sunspot (which is usually referred to as the umbra) cor- 
responds to a region of the solar surface in which the con- 
vective heat transport is impaired by an intense concen- 
tration of vertical magnetic flux. However, observations 
suggest that convection is not completely suppressed in 
these regions. In particular, there are some small bright 
features within sunspot umbrae, usually known as umbral 
dots [23l |24] , which tend to be non-uniformly distributed 
across the umbral region [25 . It is probable that umbral 
dots indicate the presence of localised magnetoconvective 
motions. These observations have motivated theorists to 
search for localised states in idealised models of magne- 
toconvection. 

It is now well-established that it is possible to find 
steady localised states (often referred to as "convec- 
tons" ) in two-dimensional Boussinesq magnetoconvection 
[26l[27]. However, the bifurcation structure of this mag- 
netoconvection system differs from that exhibited by sim- 
pler systems, such as the illustrative Swift-Hohenberg 



models, where localised states exist in a region of bista- 
bility between the conducting state and the multiple roll, 
spatially-periodic convective pattern [7l[8]. Furthermore, 
in the Swift-Hohenberg case, solution branches corre- 
sponding to localised states form a vertical "snaking" 
pattern of interlaced curves of gradually increasing num- 
bers of oscillations [28l |2^. This snaking pattern is 
mostly confined to a narrow region of parameter space, 
usually referred to as the "pinning" region. In a fi- 
nite domain, this pattern of snaking solution branches 
eventually terminates on the periodic patterned solution 
branch. In the case of magnetoconvection models, the 
presence of a conserved quantity (namely the vertical 
magnetic flux) has an important effect upon the bifurca- 
tion structure [30 . In a highly conducting fluid, convec- 
tive eddies tend to expel magnetic flux [31 , which must 
therefore accumulate elsewhere within the domain. This 
can lead to the complete suppression of convective mo- 
tions in the magnetically-dominated parts of the layer, 
leading to a flux separated state [32] in which most of 
the convective motions are restricted to field- free regions. 
This flux expulsion process enables convectons to exist in 
a subcritical regime in which the conducting state is the 
only non-localised stable solution. Numerical continua- 
tion has been used to show that the solution branches 
corresponding to localised states in magnetoconvection, 
which tend to bifurcate subcritically from the periodic 
states, exhibit a slanted snaking pattern that is not con- 
fined to a narrow pinning region [33-35 . Although it 
has been suggested that standard snaking should be re- 
covered in the limit as the box size tends to infinity 
[311 ES] 5 the existence of slanted snaking in finite domains 
illustrates the fact that the bifurcation structure of mag- 
netoconvection differs in a fundamental way from that 
of many comparable pattern-forming systems. This is 
therefore an interesting system to study in its own right. 

One of the notable features of Boussinesq magnetocon- 
vection is that it is possible to find regions of parameter 
space in which the trivial conducting state is unstable to 
a spatially-periodic oscillatory mode [22] . This raises the 
question of whether it might be possible to find localised 
oscillatory states Oscillatory convectons have been 
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found in both two [3T and three [38] spatial dimensions, 
but only in idealised models in which the vertical struc- 
ture of the system has been simplified by a projection 
onto a small number of Fourier modes. Little is known 
about the parametric dependencies of these states in two- 
dimensions and, so far, the existence of such states has 
been confirmed only in a very small region of parameter 
space (using the parameters defined below, R = 20000, 
C = 0.1, cr = 1, A = 6 and Q = 17000 - 22000 [37 ). 
Whether or not fully-resolved magnetoconvection can 
produce corresponding oscillatory states is still an open 
question. As described above, only steady localised states 
have so far been found in fully-resolved two-dimensional 
Boussinesq magnetoconvection. This also seems to be 
the case in fully-resolved three-dimensional compressible 
magnetoconvection [39^, although, again, only a limited 
exploration of parameter space has been carried out to 
date. So, on the basis of previous studies, we cannot dis- 
count the possibility that oscillatory localised states exist 
only in simplified models of magnetoconvection. Having 
said that, given the complexity of the bifurcation struc- 
ture for the steady modes it is just as plausible that these 
states do exist but have yet to be observed in numerical 
simulations. 

The primary aim of this paper is to demonstrate the 
existence of oscillatory localised states in fully-resolved 
two-dimensional magnetoconvection. We focus initially 
upon a simplified model with a truncated vertical struc- 
ture, with the aim of determining the properties of oscil- 
latory localised states in this system. We also carry out 
a brief parametric survey in order to verify that these 
oscillatory localised states are not restricted to a nar- 
row region of parameter space. So, in that sense, these 
solutions are a robust feature of this system. Having 
carried out this survey, we then move on to the more 
complex fully-resolved case. In the next section, we de- 
scribe our magnetoconvection model. Numerical results 
are presented in Section 3. Finally, our conclusions are 
given in Section 4. 



assumed to be constant, including the thermal conduc- 
tivity K:^ the kinematic viscosity z/, the magnetic perme- 
ability /io7 the thermal expansion coefficient a, and the 
(reference) density pQ. To non-dimensionalise this sys- 
tem, we follow the approach of previous authors [22 by 
scaling length and time by the layer depth, and the 
thermal relaxation time, d^//^, respectively. The veloc- 
ity field u{x^t) is therefore scaled by A^/d, whilst the 
magnetic field B (x^t) is scaled by the magnitude of the 
imposed field, \Bo\. Finally, we normalise T (x^t), which 
is a measure of the temperature of the fluid relative to 
that of the upper boundary, by the temperature differ- 
ence across the layer, AT. The dimensionless governing 
equations can then be written in the following form 

1 / Ou \ 

-(^ + {u-V)uj =-Vp + RTz 

+ (:Q{VxB)xB + V^u, (1) 
^ + {u-V)T = V'T, (2) 

-^-V x{uxB)= CV'B, (3) 

where p{x^t) is the (modified) pressure [22], whilst the 
(incompressible) velocity field and the magnetic field sat- 
isfy the constraints 

V • n = V • B = 0. (4) 

The four non-dimensional parameters in this system are 
defined in Table |l| The Chandrasekhar number, Q, is 
a measure of the dimensionless field strength. The 
Rayleigh number, i?, measures the destabilising effect of 
the superadiabatic temperature gradient. The Prandtl 
number, <j, is the ratio of the viscous to thermal diffusiv- 
ities, whilst ( gives the ratio of the magnetic to thermal 
diffusivities. It is convenient to eliminate the pressure 



II. MODEL 

II. 1. The Fully Resolved Model 

Adopting the Boussinesq approximation [22[ |40] we 
consider a layer of electrically conducting fluid that is 
heated from below, in the presence of an imposed uni- 
form vertical magnetic field. The fluid occupies a two- 
dimensional Cartesian domain of dimensions < z < d 
and < X < Ad, with periodic boundary conditions being 
enforced in the horizontal direction. The orientation of 
the Cartesian system of coordinates is chosen so that the 
z-axis points vertically upwards, which implies that the 
constant gravitational acceleration is given by ^ = —gz. 
The magnetic field is constrained to be vertical at the up- 
per and lower boundaries, which are also assumed to be 
impermeable and stress- free, and are held at fixed uni- 
form temperatures. Various properties of the fluid are 



Symbol Parameter Name Definition 



Q 


Chandrasekhar number 




R 


Rayleigh number 


ga/\Td^/ny 


a 


Prandtl number 


V 1 K 


c 


Diffusivity ratio 


r]/i<i 



TABLE I. The non-dimensional parameters in the governing 
equations for Boussinesq magnetoconvection. 



by taking the curl of Equation ([T]), which introduces the 
vorticity, u){x^t) = V x n. Given that we are working 
in two spatial dimensions, the incompressibility condi- 
tion for the velocity and the solenoidal constraint for the 
magnetic field (see Equation Q) can be satisfied auto- 
matically by introducing a streamfunction, ip{x^t): 

n = V X (V^ (x, z, t) y) = (-a^^, 0, d^^^p) , (5) 
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and a flux function, A{x^t): 

B = Bo^S/x{A{x, z, t) y) = {-d,A, 0, 1 + d^^A) , (6) 

where Bq = z is the dimensionless, uniform vertical mag- 
netic fleld HU m]. We also deflne 0{x,t), to be the 
temperature perturbation from the background conduct- 
ing state, so that 

T = 1- z^O{x,z,t). (7) 
This leads to the following set of governing equations 



dO 
dt 

dA 



d (x, z) 



d (x, z) 
d {x, z) 



■- (tV^w - aR 



00 

dx 

dV^A d [A, V^A) ' 
dz d{x,z) 



dx ' 



(8) 

(9) 
(10) 



where uj{x,z^t) is the magnitude of the vorticity vector 
(which is parallel to the ^-direction) and 



uj{x^z^t) = —V'^ip{x^z^t). 



(11) 



Note we have used Jacobian notation in equations ([8]) 
(10) where. 



d (x, z) 



d^dg_ _ d^dg_ 
dx dz dz dx ' 



(12) 



with / = f{x^z^t) and g = g{x^z^t). Having made this 
change of variables, the boundary conditions can be writ- 
ten in the following form: 



tlj = uj = = d^A = at z = 0, 1. 



(13) 



Equations ([8|-(10) have a trivial equilibrium solution 
such that 



0, 



which corresponds to a static fluid with a constant tem- 
perature gradient, permeated by a uniform, vertical mag- 
netic fleld (i.e. the basic conducting state). 



II.2. The Truncated Model 

A number of previous studies of Boussinesq magne- 
toconvection have considered a further simpliflcation to 
this model. The full horizontal structure is retained but, 
motivated by the fact that the convective pattern close to 
onset has a relatively simple vertical dependence, only a 
minimal (non-trivial) set of modes are retained in the ver- 



tical Fourier decomposition of each variable [26[ |27l |42] . 
Following these authors, we satisfy the boundary con- 
ditions (13) by considering the following Fourier expan- 
sions, 

?/^(x, t) = (x, t) sin (ttz) , (15) 
a;(x, t) = cji (x, t) sin (ttz) , (16) 
6{x, t) = Oi (x, t) sin (ttz) + O2 (x, t) sin (2tiz) , (17) 
A{x, z, t) = Aq (x, t) + Ai (x, t) cos (ttz) . (18) 

It should be noted that both the z-dependent and z- 
independent components of the magnetic flux function 
are included in order to incorporate the flux expulsion 
effect [311 [42] . It is necessary to include O2 in order to en- 
sure that the hydrodynamic problem retains some form of 
nonlinearity. Projecting onto these modes we obtain the 
following set of flve partial differential equations, which 
will become our one-dimensional truncated model [27] 

dtUJi = a (^uji — TT^cJi) — (jRO[ 

- aCQiT [(1 + A'o) {tt'Ai - A'{) + A'^'A^] , (19) 



dt02 
dtAo 



+ ^[{1+7T02)+-^10'2, 
Ifl - Ip'ldl) , 



it 



cA'^ + ^{^,A^y, 



dtAi = C (A'l - TT^Ai) + TT^i (1 + A'o) , 



with the condition that 



(20) 
(21) 
(22) 
(23) 

(24) 



(where primes denote dx)- 



U.S. Code Details 



Due to the complexity of equations (|8[)-(|1Q[) and (19) 



(14) (23) we must solve them numerically. This is accom- 



plished by discretising the equations onto one and two- 
dimensional Cartesian meshes, with grid resolution 256 
and 256 x 64 grid points respectively. Both codes are 
pseudo-spectral, using Fast Fourier Transforms (FFT) 
from standard FFTW libraries to calculate all horizontal 
derivatives. In addition, the fully-resolved model uses 
fourth order flnite differences to calculate all vertical 
derivatives. The discretised systems are evolved using a 
fourth order Runge-Kutta method for time-stepping with 
the timestep size constrained via the Courant-Friedrichs- 
Lewy (CFL) condition [47. The vorticity and the stream 
function are related by Equation ( 11 ) in the fully-resolved 



model, and by Equation ( [24| ) in the truncated model. So, 
whenever the vorticity is evolved in time, these equations 
must be inverted in order to flnd the streamfunction. For 
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the truncated model, this is accomphshed by moving into 
Fourier space, where a second derivative corresponds to 
multiphcation by (i/c)^, where /c = 27r/A is the wavenum- 
ber and A is the wavelength. Thus in Fourier space 



(25) 



where ipi and cJi correspond to the Fourier transforms 
of ipi and UJl respectively. For the fully-resolved case, 
we invert Equation (11) using an LU-decomposition [43] . 
having already simplified the horizontal dependence of 
the problem by again moving into Fourier space. 
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III. NUMERICAL RESULTS 



In this section we investigate the properties of oscilla- 
tory localised states in this magnetoconvection model, 
focusing initially upon the truncated model (in Sec- 
tion III.l) before moving on to the fully-resolved two- 



dimensional case (in Section III.2). Before proceeding, 
we first define two important quantities. In what follows, 
we shall often refer to the Nusselt number: 



N 



^0 



(1 - dO/dz) dx, 



(26) 



which measures the ratio of convective to conductive 
heat transfer [22 . We always measure this (often time- 
dependent) quantity at the base of the domain. Note that 
for N = 1 we have a purely conductive state and thus no 
convective motions. For > 1, convection is occurring, 
with larger values corresponding to more vigorous con- 
vective motions. For oscillatory states we shall see that 
the Nusselt number varies with time about a mean, given 
by TV. As described in the Introduction, convectons are 
associated with the phenomenon of fiux expulsion. To 
quantify this effect, we define 



Qeff = Q^z 5 



(27) 



where Bz is the (uniform) field in the non-convective 
part of the layer. This is a measure of the effective field 
strength in the magnetically-dominated regions within 
which the localised state is embedded. When fiux expul- 
sion is very efficient, which is certainly the case for steady 
convectons, this definition is (approximately) equivalent 
to 



geff = QAV(A-A)', 
where A represents the cell width [27^ 



(28) 



FIG. 1. (Color online) Two supercritical solutions obtained 
for the same parameter values, R = 20000, ( = 0.1, a = 1, 
A = 6 with Q = 22000, for the truncated model, indicating 
the bist ability of these states; (a) a steady convecton solu- 
tion with N = 1.31 and A = 0.996 (Qeff = 31631) (b) 12-roll 
oscillatory state with N — 1.08. Each plot shows the temper- 
ature perturbation (top) and magnetic field strength, \B\'^ 
(bottom). 



III.l. The Truncated Model 



III. 1.1. General Properties of the Oscillatory States 

Initially, we fix most of the parameters (setting R = 
20000, cr = 1, C = 0.1 and A = 6) but vary the Chan- 
drasekhar number, Q. It has been demonstrated that 
both steady and oscillatory localised states exist in the 
truncated model in this parameter regime [37], so this is 
a natural starting point for the present parameter sur- 
vey. We fix <j = 1 throughout, for numerical conve- 
nience. In a wide domain, the onset of convection in 
the absence of a magnetic field occurs at Rc ~ 657.51. 
So R = 20000 corresponds to a highly supercritical 
value of the Rayleigh number for hydrodynamic con- 
vection. Of course convection will be completely sup- 
pressed if the Chandrasekhar number exceeds some crit- 
ical value. As Q is decreased from a very large value, 
the trivial conducting state eventually becomes unstable 
at Q = Qmax = 24047 to an oscillatory convective mode 
(the initial bifurcation is always oscillatory provided that 
C < 0.83). The conducting state is also linearly-unstable 
to steady convective modes, although this occurs at much 
lower Q. Decreasing Q, the first steady bifurcation oc- 
curs at Q = QmL = 1374 <C QmL, so these bifurcations 
are not discussed further here. Henceforth, we define the 
subcritical regime to be the region of parameter space in 
which the trivial state is stable to convective perturba- 
tions (i.e. Q > Qmax), whilst the supercritical regime 
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FIG. 2. A bifurcation diagram indicating the stability range of a selection of stable steady (s) and oscillatory states (o) for 
the parameter values R = 20000, = 0.1, a = 1 and A = 6, where the dashed line (- -), at A/" = 1, indicates the unstable region 
of the trivial state. The numbers next to the branches refer to the number of convective rolls. 



is defined by Q < Qmax- Previous studies suggest that 
steady localised states exist only in large aspect ratio do- 
mains, but A = 6 appears to be sufficient for this choice of 
parameters [37 . We will investigate the effects of varying 
A in the next subsection. 

Before considering the properties of any oscillatory lo- 
calised states that might be present, we first briefly re- 
view the properties of the steady convectons that exist in 
this system. These can be located by tracking solution 
branches [26], starting from (for example) a steady mul- 
tiple roll, flux- separated state, which can be obtained by 
evolving the system from the trivial (conducting) state 
at comparatively low Q. Using a flux-separated state as 
an initial condition, the value of Q is gradually increased 
(typically using steps of Q ~ 500 — 1000), leading to 
progressively fewer convective rolls, until a single steady 
localised state is obtained. Note that it is necessary to 
take smaller step sizes (of Q ^ 100) close to the upper 
and lower stability boundaries of the convecton branches 
due to the fact that these states are more sensitive to 
small variations in Q in these regions of parameter space. 
Each state was allowed to settle for approximately 1000 
thermal diffusion times, or until the Nusselt number ex- 
hibited no further time-dependence. One of these con- 
vecton states is illustrated in Fig. [T](a). The upper plot 
in this image shows the temperature perturbation, and 
the lower plot shows the magnetic field strength, where 
lighter colours correspond to warmer fluid and stronger 
fields and darker colours correspond to cooler fluid and 
weaker fields. Steady convectons are very efficient at ex- 
pelling magnetic flux from their interior, which causes 
the effective field strength of the outside region to in- 
crease. In the solution that is illustrated in Fig. [ija), 

Qeff = 31631 > QmL. Although the convecton branch 
does extend into the subcritical regime for this particular 
choice of parameters, the solution that is illustrated here 
is supercritical and exhibits bistability with the multiple 
roll oscillatory branch, as shown in Fig. [lib). This regu- 



lar time-dependent solution takes the form of a standing 
wave oscillation in which the cells reverse their sign of 
vorticity every half period with approximately 16.2 com- 
plete oscillations per diffusion time. The reversal of the 
vorticity continuously redistributes the magnetic field, 
regularly moving the flux sheets at the top and bottom 
of the layer into regions in which the flow converges. 

Using the branch-tracking procedure, it is possible to 
locate a large number of different stable solutions in this 
system [26, 37 . Fig. |2] shows the stability boundaries for 
a selection of steady and oscillatory solution branches for 
these parameters, plotting the (average) Nusselt number, 
iV, on each branch as a function of Q. Note that we do 
not speculate on the location of the unstable solution 
branches p6 , since the location of these cannot be de- 
termined with any degree of accuracy using this partic- 
ular method. This bifurcation diagram is by no means 
complete, with numerous other solution branches exist- 
ing. These have been omitted from this plot in order 
to highlight the branches of interest. In the context of 
present study, it is worth highlighting the presence of 
several solution branches corresponding to states which 
exhibit oscillatory behaviour. These are found at lower 
values of Q than the spatially-periodic multiple roll oscil- 
latory state that is illustrated in Fig. [TJb). Indicated in 
the bifurcation diagram are 1, 2, 3 and 4-roll oscillatory 
solution branches as well as the 12 roll spatially-periodic 
state. 

Fig. [3] shows a time sequence of states (with time in- 
creasing downwards), illustrating one half-period of os- 
cillation for a 2-roll oscillatory state at Q = 18000. In 
this sequence, both rolls have the same sign of vorticity 
and oscillate in synchronisation. The sequence proceeds 
with both rolls decaying to a transitionary state in which 
the convective efficiency is lowest and the magnetic flux 
within the cells begins to redistribute. The flux sheets 
at the top and bottom of the layer then switch to op- 
posite sides of the cell and, as the convective efficiency 





FIG. 3. (Color online) A time- sequence of states showing 
a half-period of oscillation (approximately 0.059 diffusion 
times) for a two roll oscillatory solution at = 20000, cr = 1, 
= 0.1, A = 6 and Q = 18000 (Qeff ^ 25480). The snapshots 
were taken at the peaks and troughs in the Nusselt number (a) 
t = 993.396 (b) t = 993.425 (c) t = 993.455 (See Fig.[4|. Each 
of the three sets of plots shows the temperature perturbation 
(top) and magnetic field strength, \B\'^ (bottom). 



FIG. 5. (Color online) A time-sequence of states showing 
a half-period of oscillation (approximately 0.048 diffusion 
times) for a single roll oscillatory solution at = 20000, 
cr = 1, C = 0.1, A = 6 and Q = 22200 (Qeff = 25247). The 
snapshots were taken at the peaks and troughs in the Nusselt 
number (a) t = 990.002 (A^max) (b) t = 990.025 (Nmin) (c) 
t = 990.050 (A^max). Each of the three sets of plots shows the 
temperature perturbation (top) and magnetic field strength, 
(bottom). 
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FIG. 4. (Color online) A time sequence of the Nusselt number 
for the two roll oscillatory state in Fig. [3] indicating the mod- 
ulated nature of the solution. The average Nusselt number 
for this state is N — 1.25. 



FIG. 6. (Color online) A time sequence of the Nusselt number 
for one full diffusion time for the localised oscillatory state 
in Fig. |5] There are approximately 10.4 full oscillations per 



diffusion time with A^rr 



1.19, Nrr. 



1.05 and N = 1.12. 



begins to increase, they both develop into convective rolls 
with oppositely-signed vorticity to that with which they 
started. This illustrates one half period of oscillation, 
corresponding to 8.6 complete oscillations per diffusion 
time. This solution, like all of the flux separated oscilla- 
tory states, is confined to a finite range of values for Q 
(14500 <Q < 19000). 



It is also worth noting that the multiple roll flux sep- 
arated states, although they have a clear period of os- 
cillation, are modulated over much longer time-scales. 
Fig. |4] shows the variation of A/", with time, for the state 
in Fig. |3j over 15-diffusion times. We observe that there 
are distinct phases, marked by the vertical dashed lines (- 
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-), that this state transitions through, corresponding to 
different types of oscihation. Phases one {t ^ 985 — 986.9) 
and three {t ~ 992 — 994.5) correspond to the state shown 
here with the cehs oscihating in synchronisation, with the 
same sign of vorticity. Phases two {t ^ 988.8 — 991.5) and 
four {t ^ 995.8 — 999.2) correspond to a state in which 
the two cehs oscihate in synchronisation but with oppo- 
site vorticity to each other. The transitionary regions 
between these states correspond to phases in which the 
oscihation of the cehs is no longer synchronised. The 
3 and 4-ron states oscihate with similar modulated be- 
haviour. 

Taking the 2-roll oscillatory state at Q = 18000 and 
increasing the field strength, one can suppress one of the 
convective rolls so that only a single localised oscillatory 
cell remains. As indicated in Fig. |2] the 2-roll solution 
branch loses stability at Q ~ 19000. By branch track- 
ing, it can be shown that the lo branch is stable in the 
range 17900 < Q < 22200, in close agreement with the 
findings of [37 for this particular parameter regime. The 
localised oscillatory state at Q = 22200 (Qeff = 25247) is 
illustrated in Fig. [5] The cell oscillates in the same way 
as the 2-roll state, oscillating backwards and forwards, 
reversing the sign of the vorticity every half period (ap- 
proximately 0.048 diffusion times). The oscillations in 
the Nusselt number are very regular (see Fig. |6|, with 
no modulation over large time-scales. Given that the 
Lorentz force plays a crucial role in driving the oscilla- 
tions, it is instructive to look at the time-evolution of 
the spatial distribution of the Lorentz force. This is cal- 
culated as the sum of the terms proportional to Q on 
the right hand side of Equation (19). Fig. [7]is a plot of 
the spatial-dependence of the Lorentz force terms for the 
snapshots shown in Fig. [5] At each instant in time, the 
Lorentz force does appear to be nearly symmetric about 
X = 3, although there are some small asymmetries in 
the distribution (visible primarily in the narrow spikes at 
the edge of the convecton in the upper and lower plots). 
The symmetric nature of the Lorentz force distribution 
is probably due to the simplified nature of this model. 

We can see, from Fig.[7| that the oscillatory state is not 
completely flux expelled, as the Lorentz force is non-zero 
within the cell. This is in contrast to steady convectons, 
which are very efficient at expelling almost all magnetic 
ffux from their interior. Within a cell that is completely 
devoid of magnetic ffux, Aq = —1, so that all terms that 



make up the Lorentz force in Equation (19) approximate 
to zero [27 . However, for the oscillatory convecton, in- 
dicated in Fig. |5|a), at the peak of oscillation (A/'max), 
when we expect a greater proportion of the ffux to be ex- 
pelled, we ffnd that Aq ^ —0.67, indicating that the cell 
is not completely ffeld free. This is presumably due to the 
fact that convection is rather weak (A^max = 1.19). The 
inefficiency of ffux expulsion is almost certainly essential 
for the persistence of oscillations in these localised states, 
and also presumably explains their susceptibility to con- 
vective perturbations to the outside region of the cell, 
leading to their narrow range of stability. More specif- 
ically at Q ^ 22300 the state described above becomes 
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FIG. 7. (Color online) The corresponding Lorentz force plots 
for the snapshots in Fig. [5] The snapshots were taken at the 
peaks and troughs in the Nusselt number; (a) t = 990.002 

(iVmax) (b) t = 990.025 (Nmin) (c) t = 990.050 (iVmax). 



unstable, and the solution transitions to the 12-roll oscil- 
latory solution indicated in Fig. [TJb). As Q is decreased 
along the lo branch, although the motion is still domi- 
nated by a single oscillatory roll, weak convective motions 
start to appear within the magnetically-dominated region 
of the domain. This occurs at approximately Q ~ 20500 
(at which point Qeff ^ QmL). As Q is decreased further 
the state becomes increasingly perturbed and the state is 
eventually overcome by these oscillations at Q ^ 17800 
(Qeff = 19852), transitioning to the 2-roll state indicated 
in Fig. [3) Fig. [8] shows the solution right at the lower end 
of the lo solution branch. 



III. 1.2. Varying the system parameters 

In this section we assess the effects of varying the 
Rayleigh number, i?, the diffusivity ratio, and the as- 
pect ratio, A. The aim here is to determine whether or 
not these localised oscillatory states are restricted to the 
particular parameter regime that was discussed above. 
This is the ffrst time that the parameter-dependence of 
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1 2 3 4 5 6 

X 

FIG. 8. (Color online) A snapshot of a solution on the sin- 
gle roll oscillatory solution branch at = 20000, a = 1, 
C = 0.1, A = 6 and Q — 17900, clearly indicating the per- 
turbations to the outside region of the cell that develop as 
Q is decreased. This state was taken at an irregular spike 
in the Nusselt number {Nma^ = 1.34) in order to highlight 
the eddies in the outside region at the peak in their intensity 
(N = 1.20). The plot shows the temperature perturbation 
(top) and magnetic field strength, (bottom). 




X 



FIG. 9. (Color online) A snapshot of a solution on the single 
roll oscillatory solution branch at = 20000, a = 1, C = 0.1, 
A = 3 and Q = 19800 (Qeff = 22910), indicating that in 
a restricted domain the structure of the localised oscillatory 
cell remains consistent with that in a wider domain. This 
state was taken at the peak in the Nusselt number (A^max = 
1.21) where N — 1.13 for this state. The plot shows the 
temperature perturbation (top) and magnetic field strength, 
(bottom). 



these states has been explored. In fact, these solutions do 
appear to be reasonably robust, in the sense that similar 
solutions can be found across a range of different param- 
eters. In Table [IT] we summarise some of the findings 
of this parametric survey. Various properties of the os- 
cillatory state are given for Qmax (the upper end of the 
solution branch) : A corresponds to the mean width of the 
cell, whilst A'max and N give the peak and mean values of 
the Nusselt number (respectively). Also indicated is the 
location of the linear stability boundary for the trivial 
state, Q^^L M- 

One of the key findings of this survey is that the exis- 
tence of the oscillatory convecton is found to be critically 
dependent upon (. Fixing R = 20000, a = 1 and A = 6 
we find that these states can be found only for ( < 0.6. 
It is found that as ( is increased, there is a reduction in 
the stability range of these states, as indicated in Table 
[m Steady convectons exhibit a similar ( dependence [27 . 
Where oscillatory localised states do exist, the cell width 



TABLE II. A summary of the stability range of oscillatory 
convectons for a range of values of A, C and R. Indicated are 
a number of properties of the cells at Qmax: A corresponds 
to the mean width of the cell, whilst A^max and N give the 
peak and mean values of the Nusselt number (respectively). 
Also indicated is the location of the largest oscillatory Hopf 
bifurcation Qmix- Note that the period of oscillation is given 
in terms of diffusion times. 
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appears to be relatively insensitive to the precise choice 
of C, although we see a marked increase in the period of 
oscillation at higher values of (. This is presumably due 
to the fact that the lower ( oscillatory convectons are sta- 
ble at higher values of Q, which increases the amplitude 
of the Lorentz force in the boundary layers. 

For steady localised states the effect of increasing the 
box width is to shift the stability range to higher values 
of Q [37]. We observe the same phenomenon here for 
the oscillatory convectons, although the actual branch 
width is always relatively small. For steady convectons 
this branch moves completely into the subcritical regime 
for A ^ 1, with the lower end of the steady convecton 
branch positioned at Qmix |2S]- The reason for this is 
that the cell is so small in the wide domain that the in- 
crease in flux in the outside region is not large enough 
to suppress all convective perturbations if Q < Qmlx- 
In this oscillatory case, the upper limit of the stability 
range (in Q) of the localised state seems to be tending 

towards Qmax as the box size increases, with almost no 
variation in the branch width, as can be seen from Ta- 
ble |ll| Furthermore, we find no evidence for the existence 
of subcritical oscillatory convectons. This is specifically 
highlighted by the A = 16 case in Table |TI| which indi- 
cates that the solution decays rapidly to the trivial state 
almost exactly at Qmax- Unlike steady convectons, oscil- 
latory convectons can also be found in much smaller do- 
mains. As indicated in Table [HI the localised oscillatory 
state for A = 3 (which is illustrated in Fig. [9| is stable in 
the range 14400 <Q < 19800. At Q = 18000, TV = 1.25 
and A^max = 1-41, both of which compare very favourably 
to the A = 6 two roll state that is shown in Fig. [3) Sim- 
ilarly we may compare the Nusselt number for a single 
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roll oscillatory state in a domain of size A = 6/3 = 2 
and A = 6/4 = 3/2 which oscillate about average Nus- 
selt numbers of N = 1.36 and N = 1.51 respectively, in 
agreement with those for the three {N = 1.37) and four 
roll (TV = 1.50) oscillatory states in the A = 6 domain. 
We conclude from this that the 2o, 3o and 4o states in 
the A = 6 domain could be regarded simply as collections 
of interacting oscillatory convectons (any of which could 
exist in isolation). 

Varying the Rayleigh number at fixed ( (which is 
equivalent to varying the strength of the convective driv- 
ing forces) we find localised oscillatory convectons for 
5000 < R < 50000, however, for larger values of R 
(e.g. R = 70000), stable states could not be found. 
For R = 50000, oscillatory convectons are stable in the 
range 47600 < Q < 51000, however, even at the upper 
boundary (Qmax) of the solution branch {Q ^ 51000, 
A^max = 1.42, TV = 1.28, Q^^L = 68800), smah counter 
rotating eddies either side of the main cell are present, 
indicated by the fact that Qeff ^ 56848 < Qmax, similar 
to the state in Fig. [8) These perturbations outside of the 
cell for higher Rayleigh numbers are most likely the result 
of the inability of these states to expel significant quan- 
tities of magnetic fiux. As a result, significant increases 
in the thermal forcing mean that the field is not locally 
strong enough to suppress all convective perturbations 
and stabilise the layer. This could give a possible expla- 
nation as to why these states have not been observed at 
higher Rayleigh numbers. 



III.2. The Fully Resolved Model 

Having analysed the truncated model in some de- 
tail, we now turn to the model of fully-resolved two- 
dimensional magnetoconvection. In any model of this 
type, the presence of slowly decaying transients implies 
that long time integrations are required in order to ensure 
the persistence of localised states, so that the addition of 
the extra spatial structure makes it impractical to carry 
out detailed parametric surveys of the fully-resolved sys- 
tem. Therefore our parameter choices must be guided 
by the results from the truncated model. However, it 
is not always the case that all behaviour observed in 
a simplified model translates to the fully-resolved case. 
For example, consider the R = 20000, ( = 0.1, a = 1 
and A = 6 case. In the case of steady convectons, sim- 
ulations show a very dramatic increase in the stability 
range, from 13700 < Q < 24300 in the truncated model 
to 6000 < Q < 25800 in the fully-resolved model. An- 
other trait of the fully-resolved simulations is that for a 
particular set of the parameters, the cell size and convec- 
tive efficiency is significantly increased in comparison to 
the truncated model. The increased cell size can be seen 



in Fig. 10, which shows a fully-resolved steady convecton, 
found by branch tracking, as well as a 12-roll oscillatory 
state, for Q = 22000. This should be compared with 
Fig. [l] which shows the results from the truncated model 
for the same set of parameters. The Nusselt number in 
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FIG. 10. (Color online) Two supercritical solutions obtained 
for the same parameter values, R — 20000, = 0.1, a = 1, 
A = 6 with Q = 22000, for the fully-resolved model, indicat- 
ing the bistability of these states; (a) a convecton solution 
with N = 2.38 and A = 1.94 (Qes = 48048) (b) a 12-roll 
oscillatory state with N — 1.09. Each plot shows the temper- 
ature perturbation (top) and magnetic field strength, \B\^ 
(bottom). 



this case is TV = 2.38 (as opposed to = 1.31), an in- 
crease which can mostly be attributed to the increased 
cell size. Nonetheless, despite the apparent differences 
in the solutions, the presence of a steady convecton for 
these parameter values suggests that it is not unreason- 
able to suppose that the truncated model should pro- 
vide a reasonable indication of the appropriate regions 
of parameter space to search for these solutions in this 
fully-resolved model. 

Initial attempts were made to locate oscillatory con- 
vectons for R 20000, a = 1, C = 0.1 and A = 6. How- 
ever, branch-tracking by increasing Q through the ffux- 
separated oscillatory states failed to produce the required 
solution. Tracking down from higher Q (in the other di- 
rection), we found that the spatially-periodic oscillatory 
state transitions directly to a 2-roll oscillatory state. So 
we were unable to find a stable oscillatory convecton by 
following this branch-tracking process. However, it was 
noticed that there was a transient localised cell during 
the transition to the 2-roll state. Although clearly a tran- 
sient, this localised state was present for approximately 
50 diffusion times, which highlights the need for long time 
integrations in systems of this type. Given the presence 
of this long-lived transient, it is certainly possible that 
a stable oscillatory convecton solution branch does exist, 
but perhaps over a narrower region of parameter space 
than the corresponding state in the truncated model. 

Given that steady convectons tend to be wider in this 
fully-resolved model, it is not unreasonable to assume 
that the stability of these solution branches might de- 
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FIG. 11. (Color online) A time sequence of convective states 
for a two-roll flux separated oscillatory state, showing one 
half period of oscillation (approximately 0.034 diffusion times) 
at = 20000, cr = 1, C = 0.1, A = 8 and Q = 21000. 
The snapshots were taken at equally spaced time-intervals; 
(a) t = 1543.005 (b) t = 1543.022 (c) t = 1543.039. Each of 
the three plots shows the temperature perturbation, 0, (top) 
and magnetic field strength, (bottom). 
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FIG. 12. (Color online) A time sequence of states for one re- 
versal in the vorticity of the state (approximately 0.046 diffu- 
sion times) at i? = 20000, a = 1, C = 0.1 and Q = 22000. The 
width of the cell at the peak of oscillation is A = 0.787. The 
snapshots were taken at the peaks and troughs in the Nusselt 
number; (a) t = 1499.001 (A^max) (b) t = 1499.027 (Nmin) 
(c) t — 1499.047. Each of the three sets of plots shows the 
temperature perturbation (top) and magnetic field strength. 



\BY (bottom). 



pend crucially upon the aspect ratio. Increasing the box 
width to A = 8, we found a two roll oscillatory state 
at Q = 21000 (which is illustrated in Fig. [iT]). In this 
fully-resolved case, the oscillations no longer take the 
form of localised standing waves. Instead, the oscilla- 
tion has more of a horizontal component. This can be 
seen from the structure of the magnetic boundary layers 
which move gradually across the cell. The translational 
velocities are to the right in the left hand cell and to 
the left in the right hand cell. As was the case in the 
truncated model, the Nusselt number variation of this 2- 
roll state is modulated over a large time-scale, and there 
are distinct phases (of finite duration) in which the cells 
oscillate in a synchronous way. 

Taking the 2-roll state, we increased the field strength 
and found that a single roll oscillatory localised state sta- 
bilised at Q = 22000 (Qeff = 24127), as shown in Fig.fl2 



1.15 



To verify that this state was not a transient phenomenon, 
it was evolved for over 1000 diffusion times, carefully 
monitoring the mean Nusselt number, which maintains 
a stable value of N = 1.09. This average Nusselt num- 
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FIG. 13. (Color online) Time-evolution of the Nusselt num- 
ber for the localised oscillatory cell shown in Fig. |12| over 
approximately one full diffusion time. Note that there are 
approximately 12.049 full oscillations per diffusion time with 
Nme.^ = 1.13, N^in = 1-05 and N = 1.09. 
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ber is in fact slightly lower than that for the comparable 
state in the truncated model. The actual time-evolution 
of the Nusselt number is less trivial than in the corre- 
sponding solution in the truncated model. As shown in 



Fig. [13) whilst still periodic, it is clear that the Nusselt 
number has a more complex functional form. Further- 
more, the maximum Nusselt number (A/'^ax = 1.13) is 
comparatively large. This behaviour is a consequence of 
the fact that the mode of oscillation in the fully-resolved 
case does not take the form of a localised standing wave, 
which implies that there is an asymmetry between the 
vertical distribution of the hot and cold regions of the 
oscillatory convective motion. The Lorentz force, which 
is calculated as the magnitude of the third term on the 
right hand side of Equation ([T]), is indicated in Fig. 14 
for this state. The Lorentz force is calculated at the mid- 
layer of the box {z = 0.5), where the vorticity takes its 
maximum value, in agreement with the dominant modes 
of the truncation, so that the Lorentz force can be di- 
rectly compared between the two models. We observe an 
asymmetry in the profile here which is in contrast to the 
very symmetric form in the truncated model and is most 
likely responsible for the differences between their modes 
of oscillation. Due to the numerical complexity of these 
simulations we have not carried out a detailed parametric 
survey for these fully-resolved localised states. However, 
we have found that these states exist in the (approxi- 
mate) range 19000 <Q < 22000. 



IV. CONCLUSIONS 

We have demonstrated the existence of oscillatory lo- 
calised states in fully-resolved two-dimensional magneto- 
convection. This is the first time that a localised oscilla- 
tory state of this kind has been found in a fully-resolved 
magnetoconvection simulation. The time-dependence of 
these states is found to be more complex than that ob- 
served in simplified models, differing from a localised 
standing wave due to the presence of asymmetries in the 
Lorentz force. These states were located due to an under- 
standing of their stability properties that we gained by 
performing a parametric survey using a vertically trun- 
cated model. These oscillatory solutions are found to be 
restricted to the low ( regime, existing only for ( < 0.6. 
However, these states are found to be much more robust 
to variations in the Rayleigh number, with stable states 
certainly existing in the range 5000 < R < 50000, al- 
though there is no evidence for their existence at larger 
values. These oscillatory localised states expel magnetic 
flux less efficiently than their steady counterparts, so they 
are never entirely field free. This probably explains why 
they do not seem to exist in the subcritical regime. This 
point was highlighted by increasing the aspect ratio of 
the domain, a process which moved the upper stability 
boundary of these states closer to Qmax without a corre- 
sponding increase in the branch width. For a box with 
aspect ratio A = 16 it was found that the oscillatory 
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FIG. 14. (Color online) The corresponding Lorentz force plots 
for the snapshots in Fig. |12| The snapshots were taken at 
consecutive peaks and troughs in the Nusselt number, at the 
mid- layer of the box {z = 0.5); (a) t = 1499.001 (A^max) (b) 
t = 1499.027 (A^min) (c) t = 1499.047. 



suggesting that these solutions can never exhibit subcrit- 
ical behaviour. Thus we hypothesise that the existence 
of such states is a consequence of the finite geometry of 
the box. Having said that, reducing the box size appears 
to have no adverse effect upon the existence of these os- 
cillatory localised states. 

Of course, one of the main limitations of this study 
is that it is restricted to two spatial dimensions, and 
the localised solutions that have been found would corre- 
spond to roll-like structures in a three-dimensional sys- 
tem. These isolated rolls would almost certainly be un- 
stable to three-dimensional convective perturbations, so 
any localised oscillatory states in three spatial dimen- 
sions would probably have a more complicated struc- 
ture than these simple rolls. Time-dependent plume- 
like structures have been found in vertically-truncated 
three-dimensional Boussinesq magnetoconvection [38 so, 
given the findings of this paper, it is not unreasonable to 
assume that fully three-dimensional Boussinesq magne- 
toconvection should exhibit similar behaviour. Further 
work is needed to confirm whether or not this is the case. 
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It has already been established that steady convectons toconvection [39]. Work is in progress to determine the 
can be found in three-dimensional compressible magne- effects of stratification upon both steady and oscillatory 

localised states in two and three spatial dimensions. 
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